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Abstract 

The stability of functional differential equations under delayed feedback is in- 
vestigated near a Hopf bifurcation. Necessary and sufficient conditions are derived 
for the stability of the equilibrium solution using averaging theory. The results 
are used to compare delayed versus undelayed feedback, as well as discrete versus 
distributed delays. Conditions are obtained for which delayed feedback with par- 
tial state information can yield stability where undelayed feedback is ineffective. 
Furthermore, it is shown that if the feedback is stabilizing (respectively, destabiliz- 
ing), then a discrete delay is locally the most stabilizing (resp., destabilizing) one 
among delay distributions having the same mean. The result also holds globally if 
one considers delays that are symmetrically distributed about their mean. 
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1 Introduction 

We study the effect of the feedback function / on the stability of the zero solution of 
the functional differential equation 

x(t) = Lx t + sg(x t ;e) + eKf(x t ;e), (1) 

where x{t) G R n , x t G C = C([-t, 0], W 1 ), x t {6) = x(t+9) G W 1 , 9 G [-r,0], L : C -> R n 
is linear, e is a small real parameter, f,g G Cxi — ► W 1 have continuous second 
derivatives with respect to each of their arguments and satisfy /(0;e) = g(0;e) = 
for all e, and k G R denotes the feedback gain. It is assumed that the linear problem 
obtained by setting e = has a pair of complex conjugate characteristic values ±iu> / 0, 
and all other characteristic values have negative real parts. Equation ([T]) arises in 
the study of a Hopf bifurcation of an equilibrium solution, after rescaling the space 
variable x —* ex and the bifurcation parameter a — > ea; see e.g. [1]. In applications, 
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the problem is related to the feedback control of oscillations, or conversely, to the 
oscillatory instabilities arising from delayed feedback (e.g. [2,3]). 

The aim of the present paper is to obtain precise conditions under which a de- 
layed feedback action can stabilize or destabilize an equilibrium solution near a Hopf 
bifurcation, in particular when not all the system variables are available for feedback. 
Moreover, we are interested in the difference between discrete and distributed delays 
in the feedback. Taking advantage of being near a Hopf bifurcation, we use averag- 
ing theory in Section [2] to derive necessary and sufficient conditions for stability. The 
implications for delayed versus instantaneous feedback are investigated in Section [3j 
Section d] is devoted to a discussion of discrete versus distributed delays. 



2 Stability of the zero solution 

We introduce some notation. For e = 0, we write (P) as 



x(t) = Lxt 



dr)(9)x(t + 0), 



(2) 



where r\ is an n x n matrix whose components are of bounded variation on [— r, 0]. By 
assumption, (|2|) has a pair of characteristic values ±iuj ^ 0. By rescaling time it can be 
assumed that u> = 1 without loss of generality. Assume all other characteristic values 
have negative real parts. Let $ be an n x 2 matrix whose columns span the eigenspace 
of © corresponding to the characteristic value ±i. In particular, $ can be chosen such 
that 

$((9) = $(0)e je , 0e[-t,O] (3) 

where 

' -1 
1 

Similarly, let $ denote an n x 2 matrix whose columns span the eigenspace corresponding 
to H for the adjoint equation 

i-O 



J 



(4) 



drj ' (9)z(t - 6) 



(5) 



on the space C* = C([0, r], M. n ). Let F and G be n x n matrices, with elements of 
bounded variation on [— r, 0], such that 

[£>!/((); O)]0= f dF{6)m (6) 



[£>1<7(O;O)]0 

We define the scalar functions 

= tr ^ T (0)F(6')$(0) 

gx{9) = tr (V(O)G(0)$(O) 



dG(6)cj>(e)- (7) 



f 2 (9) = tr (* T (0)F(#)<D(0) Jj , (8) 
g 2 (9) = tr U T (0)G(9M0) j) , (9) 
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where "tr" denotes the matrix trace, and define the real numbers q,p by 

fO 



cos9dg\(6) + J sm6dg 2 {d) (10) 



r0 r0 

p= cos0cZ/i(0) + / sm6 df 2 (9) (11) 

Then for sufficiently small e, the stability of the zero solution of ([T]) is given by the 
following result. 

Theorem 1. Let k G M. There exists Eq > such that for e £ (0,£o) the origin is 
asymptotically stable (unstable) if 

q + Kp<0 (> 0). 

Proof. By assumption, the linear system ([2]) and the adjoint system ([5]) each have a 
2-dimensional center subspace, which are spanned by the columns of the matrices 
and ^f, respectively. The bilinear form 

(V, <p) ■■= V T (omo) - y° J° v T (C - 0)<fy(0MO <*C, (12) 

where ip G C* and 99 € C, allows a decomposition of the space C = C([— r, 0), M. n ) [4]. 
Accordingly, the solution of the perturbed equation (pQ) can be written as 

x t = <f>y(t) + X t, y(t) = (%xt) 

for some Xt G C, where y satisfies 

y(t) = Jj/(t) + e T (0) [g($y(t) + Xt! e) + + Xt; e)] , 

with J is as defined in ([J]). The change of variables y = exp(Jt)u, gives 

u{t) = ee- Jt * T (0) {g($e Jt u(t) + X f,e) + K f^ ie Jt u(t) + X f, ej) • (13) 

By averaging, one obtains the equation 

u = eg{u) +£Kf(u), (14) 



where 



/"(n)=lim^/ e- Jt y T (0)f($e Jt u;0)dt (15) 

1 /" T 

= lim - / e- Jt $> T (0)g($e Jt u;0)dt. (16) 
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It follows by the assumptions on / and g that /(0) = g(0) = 0; thus the origin is 
an equilibrium point of the system of equations (|14|) . The linear variational equation 
about the origin is 

u = E (G + kF)u, (17) 
where the averaged matrices F,G £ M 2x2 are defined by 

1 f T r° 
F= lim - / e _Jt tf T (0) / dF(6)<S>(0)e Jt dt 

G= lim ~ f e~ Jt ^ T (0) / dG{6)${6)e Jt dt, 

T^oo T J J_ T 



with F and G given in Applying Lemma 1 in [5] to F and G, we obtain 

F=^tr(^ T (0)J dF(9)$(9)^ •/- i tr ^J^ T (0) J dF(6)$(6)\ ■ J (18) 
G = ^trU> T (0)J dG(6)$(0)\ ■ I - ^tr (jV T (0) J dG(9)$(0)\ ■ J. (19) 

Prom ([3]), (jBJ), and (fTHj) . the real parts of the eigenvalues of F are both equal to 

~tr (V(o) J°dF(ewe)) 

=^tr (V(0) J° dF{e)$(0)e je ^ 
^ti(^ T (0)J dF(9)$(0)(Icos6 + Jsml 



2 
_p 
2 

Similarly, the real parts of the eigenvalues of G are equal to q/2, so that the real parts 
of the eigenvalues of the matrix G + kF in (fTTj) are given by \{q + Kp)- If (q + np) ^ 0, 
the averaging theorem implies that there exists Eq > and an almost periodic solution 
x*(e) of the original equation (JTJ) for each e 6 [0, Eq], which has the same stability 
type as the zero solution of (fTTl) [4]. Furthermore, x*(0) = 0, and x* is unique in a 
neighborhood of € C and e = 0. It follows that x*(e) = for < e < eq, since zero 
is an almost periodic solution of the averaged equation (I14p for all e. The theorem is 
then proved since the stability of the zero solution of (|17p is determined by the sign of 
q + ftp. □ 

On the basis of the above theorem, we say that the feedback is stabilizing or desta- 
bilizing depending on whether up is negative or positive, respectively. In this sense, 
the quantity p quantifies and compares the (de)stabilizing effect of the various feedback 
schemes given by /. 
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3 Delayed versus instantaneous feedback 



We now consider the role of delays in feedback with partial state information. For this 
purpose, we assume that the linearized feedback F given in (J6j) has the form 

F(9) = Ch{9) (20) 

where h : [— r, 0] — ► K is a function of bounded variation representing a scalar distri- 
bution of delays, and C £ M nxn is a structure matrix. The feedback gain k is set to 
1, or alternatively subsumed into the matrix C. The interesting case is when C does 
not have full rank, for instance when some of the system variables are not available for 
feedback. In the following we fix the number q given in ()10p by fixing L and g, and 
investigate the effect of the structure matrix C and the delay distribution h on stability. 
Let C = ^ T (0)C$(0). From PJ and (|20j), it is seen that 

P = ati{C) + f3ti{CJ), (21) 

where 

/•O [0 

a= cos6dh(9), = sm9dh{9). (22) 



In the absence of delays, i.e., when h(9) is a Heaviside step function at zero, one has 
a = 1 and /3 = 0, yielding p = tr(C) for undelayed feedback. Hence, delayed feedback 
is more stabilizing than undelayed feedback if 

atr(C) +/?tr(CJ) < tr((7) 

or equivalently, if 

(l-a)tr(C) > pti(JC) (23) 



Similarly, delayed feedback is more destabilizing if (|23p holds with the inequality re- 
versed. Although in applications the delays are often viewed as destabilizing factors, the 
condition (|23p shows the role of delays in inducing stability. In particular, if tr(C') = 
then instantaneous feedback has no effect on the stability of the zero solution. This 
case occurs, for instance, if only some of the system variables are used in the feedback. 
However, if /? tr (CJ) ^ 0, then by (I23p delayed feedback of the same variables can 
stabilize or destabilize the zero solution, depending on the values of a and (5. 

To illustrate with an example, consider the classical van der Pol oscillator with 
linear feedback control 

r° 

y + e(y 2 -l)y + y = e J [c W (t + 9) + c 2 y(t + 9)} dh(9), < e « 1. (24) 
With x = (y,y), the linear equation around the origin is 

x(t) = -Jx(t) + e ( ° Q jW) + e(° M j\(t + 9)dh(9). 
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Figure 1: Stabilization of the zero solution of the van der Pol oscillator using position 
feedback, for feedback gain c\ equal to 5(+) and 7.8(x), and a discrete feedback delay 
at r = 1. Other parameters are C2 = 0, and e = 0.1. 



giving p = ac2 — 0c\. If the feedback is instantaneous (i.e. without delays), then 
a = 1 and (3 = 0, so that p depends only on the velocity feedback C2- In this case, 
the origin cannot be stabilized if velocity information is not available for feedback. 
By contrast, if the feedback is delayed, then using only position information can yield 
stability provided (3c\ > 1, by Theorem [TJ Figure [T] shows the quenching of oscillations 
and the stabilization of the origin when C2 = and h represents a discrete delay at 



In closing this section we note that we have confined our discussion to linear feed- 
back. If nonlinear terms are added to the feedback, then it is possible to further shape 
the system's behavior in addition to changing its stability. For example, it has been 
shown that a limit cycle of desired amplitude can be created [2], or the period can be 
modified within certain limitations [3]. 

4 Distributed versus discrete delays 

An interesting question in stability investigations is how various distributions of delays 
about a given mean value affects stability. In particular, one is interested in the dif- 
ference between distributed delays having the same mean delay f = j°_ r 9 dh(9) and a 



We have 




t = 1. 
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Figure 2: The densities /i^ corresponding to uniformly distributed delays about the 
mean value f and parametrized by ji. 



discrete delay at f. For example, Ref. [6] studied the stability of the Cushing equation 
with discrete and gamma-distributed delays. In the context of a first-order system, 
it has been discussed that the stability tends to improve with increasing variance of 
the delay distribution [7], and it was conjectured that a discrete delay at f is more 
destabilizing than distributed delays having mean f [8]. A further example involving 
coupled oscillators indeed showed that increasing the variance of the delay distribution 
can enlarge the stability region in the parameter domain [9]. We shall show that these 
observations are true in a certain sense for Hopf instabilities of more general systems. 
More precisely, when the delays act towards destabilizing the system, the discrete delay 
is locally the most destabilizing one among delay distributions having the same mean 
value. On the other hand, delays can also stabilize an unstable equilibrium point, as 
seen in the previous section. In this case the discrete delay is locally the most stabilizing 
delay distribution. 

To give a systematic study on the effect of delay distributions, we let f be fixed and 
consider a family of distributions having mean value f . For this purpose, let h be some 
reference distribution with compact support satisfying^ 



dh(0) = l, / 6dh(6) = f, / {0 - tY dh{6) = <j\ (25) 

3 J — OO J —OO 

and define a family of distributions parametrized by fi > 0, 

h^e) = h(f + (6-f)/fi). (26) 
(Figure [2] depicts the rescaling (|26|) for the case of uniformly distributed delays.) It is 



1 For simplicity, here we view h as a probability distribution, that is, a monotone function satisfying 
(|25[) . The particular choice of h will not be important for the following discussion. 
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easy to check that 



"OO 



-oo 



/CO l-OO 
0<%(0)=r, / (8-f) 2 dh fl (9)= f i 2 a 2 
-OO J —OO 



Hence, the family hn provides a natural way to change the variance fj?a 2 of the distri- 
bution while keeping the mean value fixed. We denote the corresponding values in (I22p 
as 

/OO /"OO 
cose dh^e), ^= smOdh^e). (27) 
-00 j —00 

We also define 

«o = cos t, Pq = sinf, (28) 

i.e., the values of a, (3 for a discrete delay at r, which is equivalent to letting ho be a 
Heaviside step function at f. The question is, for a fixed structure matrix C, how the 
quantity 

Pll = a^tT(C) + frMJC) (29) 

changes as /i is varied. 

Before considering the general case, it is instructive to first look at the specific 
example of uniformly distributed delays shown in Figure El Here one has 

ct u = — sin ijl cos t, j3 u = — sin u sin r, 
A* A* 

which yields 

^ _ sin^ ^ tr ^ cog f + tr (j(7) sin -rj 

sin A* / Qn N 

= Po- (30) 

A 4 

It is thus seen that the dependence of on the parameter fi is not monotone. In fact, 
by Theorem [TJ the sign changes of for varying values of /i indicates that stability 
switches can occur as the variance of the delay distribution is changed. Nevertheless, 
since | sin/i| < for all \i > 0, one has \po\ > \pp\, which shows that a discrete delay 
has a stronger effect on stability than all uniformly distributed delays having the same 
mean value. Moreover, there are certain values of the distribution variance (given by 
sin ijl = 0, [i > 0) for which the feedback has no effect on stability. This last property is 
purely an effect of the delay variance and is independent of the choice of the structure 
matrix C. 

The extremal property of discrete delays observed above can be extended to more 
general delay distributions. We first give a local characterization. 



Proposition 2. 



dp^ 



d 2 v 

= 0, and -—4r 

M =0 d ^ 



-a 2 p 
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Proof. By a change of variables in (p7|) and using ([26]) . we have 

/OO /'OO 
cos 9 dh fM (0) = I cos(f + /x(s — f)) dh(s). (31) 
-OO J —oo 



Using (I28p it is seen that and (3^ are smooth functions of /i on 1. Differentiating 
under the integral gives 



/OO 
sin(f + p(s — f))(s — f) dh{s) 
-OO 



Thus da^/d ^1^=0 = 0. Similarly, dp^/dpl^o = 0. On the other hand, 



d 2 



a,, 



dp? 

d 2 P» 



dp 2 



-a 2 cost = — o~ 2 ao, 



-a 2 sinr = — a 2 Po- 



Using in j29J), we obtain the conclusion. □ 

By the above result, if po is nonzero, then it is a local extremum for p^. This shows 
that discrete delays are indeed special in a certain sense. Thus, if the delayed feedback is 
a destabilizing one (p^ > 0), then a discrete delay is locally the most destabilizing delay 
distribution, and increasing the variance of the distribution reduces p^. This confirms 
the observations of [7-9] and shows that it is generally true near a Hopf instability. 
However, as noted above, delays can also have a stabilizing effect (pn < 0), in which 
case a discrete delay is locally the most stabilizing one, and increasing the variance 
of the distribution can yield instability. In both cases, increasing the variance of the 
distribution locally about a discrete delay reduces the effect of delays in the feedback. 

For symmetrically distributed delays, we can also give a global characterization of 
the extremal property of discrete delays. 



Proposition 3. For delay distributions that are symmetrically distributed about their 
mean valu^, \p^\ < \po\ for all p > 0. 

Proof. Expanding the cosine term in (|3ip . 

/oo roo 
cos(//(s — f)) dh(s) — sin f / sin(/i(s — f)) dh{s). 
-oo J — oo 

The second integral vanishes because the distribution is symmetric about r and sine 
is an odd function. Thus by (|28l) . = «o cos(/i(s — f))dh(s), and similarly 
= Po cos(/u(s — f )) dh{s). Hence, from ([29]) ). 



Pfi=Po cos(/x(s - r))dh(s). 



2 That is, h'(r + 8) = h! (t — 0), where the derivative exists a. e. by assumption. 
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Then the estimate 



/OO /"OO 
I cos(^(s - f))| dh(s) <\po\ / dh(s) = \p \ 
-oo J — oo 

follows. □ 

Finally we note that po depends only on the mean delay and not on the distribution 
h. Hence, the extremal properties of discrete delays given in Propositions [2] and [3] are 
independent of the particular choice of the reference distribution h. 
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